Semiclassical quantization of the Bogoliubov spectrum 
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We analyze the Bogoliubov spectrum of the 3-sites Bose-Hubbard model with finite number of 
Bose particles by using a semiclassical approach. The Bogoliubov spectrum is shown to be associated 
with the low-energy regular component of the classical Hubbard model. We identify the full set of 
the integrals of motions of this regular component and, quantizing them, obtain the energy levels of 
the quantum system. The critical values of the energy, above which the regular Bogoliubov spectrum 
evolves into a chaotic spectrum, is indicated as well. 
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The Bose-Hubbard model (BH- model), 



i=i 



ai 



U 



L 

i=i 



(1) 



constitutes one of the fundamental Hamiltonians in the 
condensed matter theory. The number of phenomena, 
discussed in the frame of the BH-model, is so diverse that 
sometimes it is difficult to see any link between them. In 
particular, this concerns the phenomena of superfluidity 
and Quantum Chaos. Indeed, the former phenomenon 
assumes the regular phonon-like excitation spectrum, de- 
scribed by the Bogoliubov theory [l], [|J , while the latter 
phenomenon implies a highly irregular excitation spec- 
trum, described by the random matrix theory [1, 0, 
This seeming contradiction is resolved by noting that 
these two spectra refer to different characteristic energies 
of the system. It is the aim of the present work to un- 
derstand of how the regular Bogoliubov spectrum evolves 
into an irregular one as the system energy is increased. 

To approach the outlined problem we consider the sim- 
plest non-trivial case of the 3-sites BH-model. Indeed, as 
shown below, the 3-sites BH-model is well entitled for the 
Bogoliubov analysis and, at the same time, is known to be 
chaotic. It is worth of noting that the 3-sites BH-system 
has been intensively studied during the last decade with 
respect to phenomenon of self-trapping in the system of 
coupled nonlinear equations @, H[ , generalization of the 
dynamical regimes of the celebrated 2-sites BH-system 
ld |. and as a model for multiparticle quantum chaos 
In the present work we use the 3-sites BH-system 
as a model for studying the Bogoliubov spectrum of the 
interacting Bose particles in a lattice. 

First we discuss the Bogoliubov spectrum of this 
system, following the particularly suited for our aims 
method of Ref. [g]. The starting point of the analysis 
is the Hamiltonian of the BH-model in the Bloch basis, 
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which ones obtains from (fTJ) by using the canonical trans- 
formation bk — (1/VL) 'J2 l exp(z27rfcZ/L)a;. The Hilbcrt 
space of @ is spanned by the quasimomentum Fock 
states \ri—-i, no, where ^ fc rik = N is the total num- 
ber of particles. To find the Bogoliubov states of the 
system ^ one uses the anzatz Q 
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Substituting (J3j) in |2]| and assuming the limit N — > oo, 
U — ► 0, g = UN/L = const, we obtain the following 
equation on the coefficients c n , 



2(6 + g)nc n + gnc n -i + g(n + l)c n+1 = Ec n 



(4) 



where g — NU/3 and S = [1 — cos(27r/3)] = 1.5 are 
the so-called macroscopic interaction constant and single- 
particle excitation energy, respectively. (From now on we 
set the hopping matrix element T to unity.) The matrix 
equation ([4]) can be solved analytically by mapping it to a 
differential equation. Namely, introducing the generating 
function $(0) = (1/V2tt) J2 n c n e i2ne , we have 



Heff = (& + g)h + g( 
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h = -id/dO .(6) 
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It takes only a few lines to prove that the spectrum of 
the eigenvalue problem §5$ is given by E n = —E + 2Qn, 
where SI = \J 2gS + 6 2 is the Bogoliubov frequency and 
Eq = Sl/2 — (6 + g) is the Bogoliubov correction to the 
Gross-Pitaevskii ground energy Q. 

It should be mentioned that, instead of ([3]), one can 
use a slightly different ansatz, \^>) = J2 n c n\ n , N — 2n — 
p,n + p). Repeating the above described procedure, we 
again obtain a linear spectrum which, however, is now 
shifted by p quanta of the Bogoliubov energy with re- 
spect to the previous spectrum. As a consequence, the 
excited energy levels of the Bogoliubov spectrum are mul- 
tiple degenerate. To fix the notation we shall call the 
number m (which counts the Bogoliubov levels begging 
from the ground state m = 0) by the primary quantum 
number and the number j, which labels the degenerate 
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sublevels, by the secondary quantum number. This sec- 
ond quantum number takes integer values for even to and 
half-integer values for odd m in the interval \j\ < m/2. 
Thus we have 

E mJ =Qm, fl= v / 2g6 + S 2 , \j\ < m/2 . (7) 

The corresponding to ([7]) wave functions are denoted by 
\^m.j) where, by construction, \^ m ,j) — \m/2 — j, N — 
to, to/2 + j) for g = 0. 

We would like to stress that the above result is valid 
only in the limit N = oo and for any finite N the low- 
energy spectrum of the system deviates from ([7]). As an 
illustration to this statement, Fig. [JJ shows the energy 
spectrum of the 3-sites BH-model for N — 40 and < 
g < 4. It is seen in the figure that (i) the spectrum is 
not linear, and (ii) the Bogoliubov levels are splitted with 
respect to the secondary quantum number. In the rest 
of the paper we shall quantify both of these effects by 
using a semiclassical approach. This approach will also 
allow us to indicate the critical value of the energy above 
which the spectrum of the BH-model is chaotic. 

As an intermediate step, let us derive the Bogoliubov 
spectrum ([7J by using semiclassical arguments. The clas- 
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FIG. 1: Low-energy spectrum of the 3-sites BH-model for N = 
40. The energy is measured with respect to the ground energy 
Eq and scaled with respect to the Bogoliubov frequency Q. 
Error bars indicate the energy intervals, where the classical 
counterpart of the system shows chaotic dynamics. 



sical counterpart of the Hamiltonian {5J is obtained by 
scaling it with respect to the total number of particles, 
H/N — > 77, and identifying the operators b\./y/N and 
bk/V~N with pairs of canonically conjugated variables 
(6 1 , fefc) . Next we switch to the action-angle variable, 
bk — V7t ex P(*0fc)j an d explicitly take into account that 
Ik = 1 is an integral of motion. This reduces our sys- 
tem of 3 degrees of freedom to a system of two degrees 
of freedom, 

H = (5 + g)(I- 1 +I +1 ) + 2gIo^rj^cos{8- 1 +6 +1 ) 
+ 2.g£± h^/hTTi cos(20 Tl - 9 ±1 ) , 

(8) 

where Iq = 1 — J_i — 7+i and the phases 9±i of variables 
b±i(t) are measured with respect to the phase of bo(t). 
The low-energy dynamics of the system ([5]). which is as- 
sociated with the low-energy spectrum of the system (JlJ , 
implies I±\ <^ Iq. Keeping in the Hamiltonian Q only 
the terms linear on I±i, and using one more canonical 
transformation, 

/ = / + i+J_i, 9 = (0+i + 0_i)/2 , (9) 
J = (7+i - 7_i)/2 , = 0+i-0_i, 

we obtain 

H eff = (S + g)I + g\/l 2 - 4J 2 cos(26>) . (10) 

[Note that for J = this Hamiltonian coincides with the 
classical counterpart (n/N — > I) of the effective Hamil- 
tonian ([6|).] Finally, we integrate the system (fTO]) by in- 
troducing a new action, 

i = ^£ i(e,E)de , (ii) 

and resolving Eq. (fTTj) with respect to the energy. This 
gives 

E = fll, fl = ^2gS + 5 2 , 7<1. (12) 

Note that the energy is independent of the action J, 
which may be chosen arbitrary in the interval \J\ < 1/2. 
Referring to the original quantum problem, this action J 
obviously labels the degenerate sublevels of the excited 
Bogoliubov states [ll| . 

Now we are prepared to discuss the finite-iV correc- 
tions to the Bogoliubov spectrum. Let us analyze the 
dynamics of the classical system (|8|) in more detail, with- 
out assuming I±i <C 7o- First we shall consider the sym- 
metric solutions, where b—\(t) — 6+i(t). Expressing the 
Hamiltonian ((SJ) in terms of the canonical variables © 
and setting there d = and J = 0, we have 

H 1D = {5 + g)I + g{l - 7)7 cos(26») (13) 
-3gl 2 /<i + gl ^2(1 -7)7cos(0) . 

The phase portrait of the ID system (fT3"|) is depicted in 
Fig. [21a) . Our particular interest in this phase portrait 
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FIG. 2: Left: phase portrait of the ID system (|13[) for g = 1. 
The bold lines restrict the chaotic region, where the trajec- 
tories are unstable with respect to variation of (J, Right: 
slow dynamics of the variables J = (J) and $ = for g = 1 
and / = 0.133. 

are the trajectories near the origin 7 = 0, which can be 
associated with the Bogoliubov states. It is seen in the 
figure, that these trajectories are strongly affected by the 
elliptic point in the upper part of the phase space. As a 
consequence, the eigen-frequency of the system depends 
on the action I. Namely, ft — Cl(I) vanishes for the 
separatrix action I* , and for I -C I* one has 

n(i) = n-ji, (w) 

where the nonlinearity 7 is a unique function of g. (For 
instance, for g — 0.1,1,4 we have 7/fi = 0.1,0.6,1, 
respectively.) Referring to the original quantum prob- 
lem the result (fP4l means that the energy difference be- 
tween (m -I- l)-th and m-th Bogoliubov levels decreases 
as jm/N. 

Next we address the 'stability' of the symmetry plane 
trajectories depicted in Fig.fSfa) with respect to variation 
of J. Within the Bogoliubov approximation the action 
J is an integral of motion and may be chosen arbitrary 
in the interval | J| < 1/2. It should be understood, how- 
ever, that in reality the action J does depend on time. 
An example of this dependence is given in Fig. [3J It is 
seen that the time evolution of the system is a super- 
position of fast dynamics, where J(t) and $(t) oscillate 
with the Bogoliubov frequency (more precisely, with the 
frequency f2), and slow dynamics, with the characteristic 
frequency of the orders of magnitude smaller than the Bo- 
goliubov frequency. Going ahead, we note that this new 
frequency defines the splitting of the Bogoliubov levels 
in Fig. [1] and for the moment we only stress that the 
system dynamics remains regular. This conclusion holds 
for any trajectory of the effective ID system (TPS")) , pro- 
viding that the trajectory lies well below the separatrix. 
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FIG. 3: An example of time evolution of the conjugated vari- 
ables J(t) and i9(£) in the regular regime. Parameters corre- 
spond to the second from top trajectory in Fig. [2|a). 

If we choose a trajectory closer to the separatrix, we ob- 
serve a transition from regular to chaotic dynamics. We 
identify the exact border of the transition to chaos by 
calculating the Poincare cross-section of ([5]) for different 
values of the energy E and evaluating the volume of the 
chaotic component as the function of energy [12]. It is 
found that the chaotic region is restricted to a relatively 
narrow energy interval E min (g) < E < E max (g). For 
g = 1 the phase trajectories of the system (fT3|h corre- 
sponding to E m in and E max , are marked by the bold 
lines in Fig. ^a). Additionally, the error bars in Fig. [T] 
indicate the chaotic energy intervals for different values 
of g. The depicted borders are consistent with the vi- 
sual analysis of the spectrum and suggest the following 
simple criteria of the transition to chaos: it takes place 
when the total splitting of the Bogoliubov levels with re- 
spect to the second quantum number exceeds the mean 
distance ft between the levels. 

The question on the sublevels splitting is in turn. As 
it was already mentioned, this splitting is defined by the 
slow dynamics of the system. To address this slow dy- 
namics we introduce the new variables J = (l/2ir) <f JdO 

and d = (l/27r) <f fldO, where 9 = fit is the phase con- 
jugated to the action I. (Note that the action / is an 
adiabatic integral of motion and, hence, does not depend 
on time.) The Hamiltonian equations of the motion for 
the variables J and d read, 

d = (dH/M) w -2gJ, (15) 
J = -(dH/dJ) 6gV{g, I) sin(3?9/2) , 

where V(g, I) w ((//2) 3 / 2 cos#) and (. . .) means time av- 
erage over one period of the fast dynamics. Thus the slow 
dynamics is defined by the pendulum-like Hamiltonian, 

H s i ow = -gJ 2 + AgV(g, I) cos(3t9/2) . (16) 
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It is worth of noting that, to obtain (fT6]l . we have as- 
sumed the quantity V(g, I) to be independent of J, which 
can be justified only if J <C 1/2. Nevertheless, the Hamil- 
tonian (| 1 6[) is found to well capture the main features 
of the low-energy regular dynamics for arbitrary J. In 
particular, it correctly predicts the existence of stable 
points at i? = 0, ±47r/3, where the phases of b±i(t) are 
locked to and 120 degrees with respect to each other [see 
Fig. [2jb)]. The size of the stability islands around these 
fixed points is obviously given by the separatrix trajec- 
tory of the pendulum, i.e., is proportional to \V(g, I)] 1 ^ 2 - 
It is instructive to consider the limiting case g — ► 0. 
As easy to show, in this limit V(g, I) — > and, hence, 
Hsiow = —gJ 2 - Let us prove that this results cor- 
responds to the first order quantum perturbation the- 
ory on U. Indeed, calculating the first order correc- 
tions to the energies of the quasimomentum Fork states 
|* mj -) = |m/2 — j,N — m,m/2 + j), we have AE = 

(U/2L)(9 mJ \ £ blhl^hJih + k 2 - k 3 - fcOl* m j> ~ 

-(U/L)j 2 , or AE = -gJ 2 , where J = j/N and E = 
E/N. Thus for small g the splitting between sublevels 
grows linearly with g. This linear regime changes to a 
nonlinear one as soon as the second term in the Hamil- 
tonian (fl"6)) takes a non-negligible value. This second 
term also causes the rearrangement of the sublevels, 
clearly seen in Fig. [1] Needless to say that in this case 
the second quantum number is defined by the action 
J = (l/27r) § Jdd, which amounts to the phase volume 
encircled by the trajectories in Fig. H£b). 

We conclude the paper by formulating quantitative cri- 
teria for the onset of Quantum Chaos. As mentioned 
earlier, the transition to irregular spectrum occurs when 
the total splitting of the m-th Bogoliubov level compares 
with the Bogoliubov frequency. Ignoring the nonlinear 
corrections, one has g(m + 1) 2 /AN ~ il, or 

m cr « - . (17) 

V 9 J 

Through the relation E cr « f2m cr this estimate defines 
the critical value of energy above which the regular spec- 
trum transforms into a chaotic one. 

In conclusion, we have analyzed the Bogoliubov spec- 
trum of the BH-model. This spectrum corresponds to 
low-energy excitations of the system and is usually intro- 
duced by using the Bogoliubov-de Gennes transforma- 



tion. This standard method, however, is rather formal 
and hides the underlying classical dynamics of the BH- 
model. In this work we use a semiclassical method which, 
by definition, explicitly refers to the classical dynamics 
and provides in this way a deeper insight in the structure 
of the low-energy spectrum of the system. In the present 
work we restricted ourselves by considering the 3-sitcs 
BH-model, although many of the reported results hold 
for L > 3 as well. An advantage of the 3-sites model 
is that, thanks to a relative low dimensionality of the 
system, its classical dynamics can be understood in ev- 
ery detail. In particular, the phase space of the system 
essentially consists of two regular and one chaotic com- 
ponent in between, where the low-energy regular com- 
ponent is shown to be associated with the Bogoliubov 
spectrum. We identify the full set of the integrals of mo- 
tion for this low-energy regular component and, quantiz- 
ing them, obtain the low-energy levels of the quantum 
BH-model. These levels are labelled by two quantum 
numbers, to and j. The first quantum number to corre- 
sponds to usual Bogoliubov ladder, where the distance 
between neighboring levels is approximately given by the 
Bogoliubov frequency Q (i.e., E m+ \^ — E m j ~ Q). The 
second quantum number j labels (to + 1) sublevels of 
the m-th Bogoliubov level, where the splitting between 
the sublevels is proportional to the interaction constant 
g and inverse proportional to the system size N (i.e., 
E m ,j+i — E m j ~ g/N). If we go up the energy axis, the 
total splitting of the Bogoliubov levels compares the dis- 
tances between the levels and the energy spectrum shows 
a transition from a regular to irregular (chaotic) one. 

The described scenario of evolution of the Bogoliubov 
spectrum into a regular, Bogoliubov-like spectrum and 
further into a chaotic spectrum also holds for the BH- 
model with L > 3 sites. However, to indicate the critical 
energies for these transitions remains an open problem. 
The qualitative difference between the 3-sites and, for 
example, 5-sites BH-models is that the latter system has 
two different Bogoliubov frequencies, associated with two 
different single-particle excitation energies. For some val- 
ues of the macroscopic interaction constant <?, these fre- 
quencies become commensurable, which strongly affects 
the onset of chaos. We reserve this problem of interacting 
Bogoliubov spectra for future studies. 

Fruitful discussions with S. Tomsovic and partial sup- 
port by Deutsche Forshungsgemeinschaft (within the 
SPP1116 program) is gratefully acknowledged. 
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